Non-contact rack and pinion powered by the lateral Casimir force 
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The lateral Casimir force is employed to propose a design for a potentially wear-proof rack and 
pinion with no contact, which can be miniaturized to nano-scale. The robustness of the design is 
studied by exploring the relation between the pinion velocity and the rack velocity in the different 
domains of the parameter space. The effects of friction and added external load are also examined. 
It is shown that the device can hold up extremely high velocities, unlike what the general perception 
of the Casimir force as a weak interaction might suggest. 
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With the emergence of the new generation of minia- 
turized mechanical devices such as micro- and nano- 
electromechanical systems (NEMS), we have been wit- 
nessing a paradigm change in the technical problems in- 
volved in making machines and the strategies needed to 
resolve them. In particular, tribological interactions, i.e. 
friction, adhesion, and wear, appear to pose new chal- 
lenges at small length scales The abundance of sur- 
faces in contact in small devices and high friction are 
problematic, and the traditional use of lubricants does 
not work because they become excessively viscous when 
made into molecular layers Q. Devices with sliding 
surfaces in contact are known to wear out too rapidly 
Q , and it seems that strategies to minimize contact be- 
tween the surfaces are needed to help make them more 
durable. The presence of short-ranged attractive disper- 
sion or Casimir forces can cause tiny elements in small 
devices to stick together and bring the devices to stop 
4], and a line of ongoing active research is focused on 
devising novel techniques to avoid such effects @, Q . In 
light of these technical difficulties, it seems desirable to 
have novel designs for mechanical devices that can oper- 
ate without physical contact between their parts. 

As a key interaction at nano-scale, Casimir force can 
be harnessed and used in small devices, as demonstrated 
by Capasso and collaborators who developed an actua- 
tor powered by the normal Casimir force between a flat 
plate and a sphere @, 0] ■ To avoid the limited applica- 
bility that the parallel-plate geometry might offer, one 
can make the two surfaces corrugated and take advan- 
tage of a lateral component to the Casimir force, as has 
been recently proposed Q and indeed verified experimen- 
tally [10(. The lateral Casimir force between corrugated 
surfaces, provides a possibility for friction-less transduc- 
tion of lateral forces in nano-mechanical devices without 
any physical contact between them. The coupling be- 
tween two surfaces via the quantum vacuum is realized 
by a term proportional to the sinus of the phase differ- 
ence between them and is a macroscopic manifestation of 
quantum coherence, which is reminiscent of the Joseph- 
son coupling between superconductors [ij EH . Because 



the coupling is nonlinear, its mechanical response could 
involve oscillatory and unstable behavior similar to the 
what is observed in Ref. § for the case a device powered 
by the normal Casimir force. 

Here we make use of the lateral Casimir force to de- 
sign a nano-scale rack and pinion without intermeshing, 
as shown in Fig. QJi. Our system consists of a corrugated 
plate (rack) and a corrugated cylinder (pinion) that are 
kept at a distance from each other. The pinion could be 
subject to external load (see Fig. |Tb), and could experi- 
ence friction when rotating around its axis. For uniform 
motion of the rack with a velocity Vr, we find that the 
"contact" pinion velocity Vp (see Fig. QJ,) is locked in 
to the rack velocity for sufficiently small values of Vr, 
and that there is a threshold rack velocity at which the 
pinion undergoes a skipping transition where the pinion 
can no longer hold the cogs in perfect registry with the 
rack. In the skipping regime, we find that the average 
pinion velocity can be both positive and negative (with 
Vr > 0) depending on the initial phase mismatch be- 
tween the corrugations. The effects of the external load 
and friction are considered and the regions in the param- 
eter space where the pinion can do work against the load 
are determined. These studies could help us examine the 
feasibility and efficiency of such a design as a mechanical 





FIG. 1: (a) The rack and pinion with no contact, (b) The 
schematics of the rack and pinion with an external load W. 
(c) Two corrugated plates with lateral shift x — y. The equi- 
librium position is at x = y. 
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transducer. 

When two sinusoidally modulated plates with identical 
wavelengths of A are shifted with respect to each other by 
a length x—y (see Fig. QJ;), the lower plate experience the 
lateral Casimir force ^lateral = — Fsin [^-(x — y)] 0, [U 
where the amplitude F depends on the mean separation 
of the plates and the amplitude of corrugations 
13] . The lateral Casimir force introduces a net torque on 
the pinion, which plays the central role in the equation 
of motion for the coordinate x = R9 (Fig. [T|d). The 
equation of motion reads 
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where / is the moment of inertia of the pinion about its 
major axis, £ is the rotational friction coefficient, and W 
is an external load against which the pinion should do 
work. 

We focus on the uniform motion of the rack with 
a velocity Vr, i.e. y = Vnt, although similar analy- 
ses can be performed for other types of motion such 
as vibrating 14j or undulating 15] racks, in both of 
which cases the motion can be rectified. The non- 
linear equation [T] can be better studied in the phase 
plane (u = 2n(x — y)/X,v — ii), where it reads u = v, 
v = — sin it — e(v — vq) — W/F. In the above equations, we 
have measured the time in units of T = -J IX/(2itFR 2 ), 
and have defined the dimensionless parameter e = T^/J, 
which is a measure of the relative importance of friction 
in the system. The parameter vq = uq = 27r(±o — VrT)/A 
is the initial value for v. A key velocity scale is given by 
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which corresponds to the velocity at which the kinetic en- 
ergy of the rotating pinion is of the order of the work done 
by the lateral Casimir force upon displacement by one 
tooth. As we will see below, the quantity FX can be con- 
sidered as the effective "bond strength" of the coupling 
between the pinion and rack, so that a "bound state" 
between them can only tolerate pinion kinetic energies of 
this order or magnitude, which means that Eq. gives 
the skipping velocity. For simplicity, we only consider the 
case of pinions that are initially at rest throughout this 
paper (±o = 0), which means vq = — Vr/Vs, although 
the formulation can be readily used to study other initial 
conditions as well. 

No dissipation. In the absence of dissipation and 
load, Eq. Q] is identical to the celebrated nonlinear 
plane pendulum problem (ljj . It is well known that 
for this system, the "energy" h = \v 2 + 1 — cos it = 

i (jzfj + 1 — cosuo, is a constant of motion, and there 
are two families of periodic orbits, corresponding to ro- 
tations (h > 2) and oscillations (h < 2). For < h < 2, 
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FIG. 2: (a) Pinion velocity versus rack velocity in the absence 
of dissipation and external load, (b) Domains of positive and 
negative pinion velocities as well as the phase boundary for 
the skipping transition. Only half of the plot is shown due to 
mo — > — tto symmetry in this case, (c) Phase boundary for the 
skipping transition, in the absence of dissipation, as a function 
of the external load, (d) Phase boundary for the skipping 
transition, in the presence of wea k dissipation, as a function 
of the external load for \/Fu/F = 0.05. (e) Pinion velocity 
versus rack velocity in the presence of weak dissipation and 
external load, (f) Force-velocity response of the rack and 
pinion could be of two general forms depending on the value 
of the rack velocity. 



the system is oscillatory with a period of AK{^/h/2), 

where K(m) = J^ 2 d0(l — m 2 sin 2 (?) -1 / 2 is the com- 
plete elliptic integral of the first kind [13] ■ In this case 
|u| | cos _1 (l — h)\ ^ 7r, which means that the distance 
| a; — 7/| between the teeth of the rack and pinion does 
not exceed A/2. Thus the pinion is locked- in with the 
rack and will have a forward motion with superimposed 
oscillations (jerks). In the case of h > 2, the system 
is not locked in any more and depending on the initial 
conditions it can have different behaviors. 

The value of h, which determines the behavior of 
the system, depends on Vr and uq. For Vr < 
Vs \/2(l + cosito) the rack and pinion are geared up and 
we have Vp = Vr for the jerk-averaged pinion velocity. 
Higher rack velocities Vr > V~s\/2(1 + cosu ) cause the 
system to go to the "rotation" phase of the equivalent 
pendulum problem, which means that the pinion skips 
teeth with respect to the rack, but will still have a net 
average velocity that can be calculated as 
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In Fig. [2£i, the pinion velocity is plotted as a function of 
the rack velocity for two values of the initial phase mis- 
match. The pinion velocity rises linearly with the rack 
velocity initially and then drops abruptly at the skipping 
transition. At large rack velocities Vr 3> Vs, the pinion 
velocity has an asymptotic form Vp = cosuoVg 2 /VR + - ' ' > 
which shows that it vanishes at infinity, and that the de- 
cay can be both from below (for uq > 7r/2) and above (for 
uo < 7r/2). An intriguing feature here is the possibility 
of dropping into negative pinion velocities — the reverse 
gear — after the skipping transition, as Fig. [5^ shows. 
Figure [2}d delineates the different domains of phase lock- 
in, forward motion, and reverse motion, in the space of 
the rack velocity and the phase mismatch. One can see 
that for < uo < tt/2 the average velocity is always pos- 
itive, while for ir/2 < uq < ir the pinion velocity can be 
switched from positive to negative by increasing Vr/Vs- 
This can be achieved by either speeding up the rack or 
decreasing the skipping velocity by increasing the sepa- 
ration between the rack and pinion (see below). 

The effect of external load. In the absence of dissipa- 
tion, Eq. Q] has ti = v 2 /2 + 1 — cos u + Wu / F as the 
constant of motion. Two classes of motion are separated 
by the saddle-loop of energy h' s = 1 + y/l - (W/F) 2 — 
[n- sin' 1 (W/F )}(W/F), leading to rotations for ti > h' s 
and oscillations for ti < h' s , in the equivalent pendulum 
problem. Figure Gfc shows the domains for two different 
regimes as well as the boundary for the skipping tran- 
sition as a function of the rack velocity and the initial 
phase mismatch. As the load is increased, the two fixed 
point approach each other thereby narrowing down the 
locked-in region, until at W — F it disappears when they 
fully merge. Therefore, W < F is a necessary condition 
to have positive rack velocities. 

Weak dissipation. The dynamical system described 
by Eq. [1] is not integrable in the presence of the dissipa- 
tion term, as energy is not conserved. If the dissipation 
is weak, such that (Vr/R 2 + W < F, we can use the 
Melnikov method 1|| 3| to study the perturbed phase 
portrait of the system. In this regime, the force scale 
Fpj = 2 niFi 2 seems to play an important role. (Note that 

Similar to the dissipation-free case, the system devel- 
ops the two classes of commensurate and incommensu- 
rate motion transduction separated by a skipping transi- 
tion boundary. Figure [2ji shows an example of this phase 
diagram that has been calculated numerically, for differ- 
ent values of the external load. The phase boundary con- 
sists of horizontal (uo-independent) parts that are met by 
a bell-shaped central part. In the part that does not de- 
pend on u , we find that for (V R /R 2 < ^^FFh-W the 
system is in the locked-in phase and we have Vp = Vr. 
As the rack velocity is increased above this limit, the 
system undergoes a skipping transition. The pinion ve- 
locity for (Vr/R 2 > ^y/FFp, - W can be found from 
Eq. [3J with a value h — h m that is a solution to this 



equation Vr 
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E(m) = J^ 2 dOyfl — m 2 sm 2 9 being the complete el- 



liptic integral of the second kind 17[ ■ Figure ^ shows 
the pinion velocity as a function of the rack velocity for 
different values of the external load. The general fea- 
ture of a drastic drop in the pinion velocity after the 
skipping transition is observed, and the asymptotic be- 
havior V P = -WR 2 /C + k V s/ V i + ■ ■ ■ , at large rack 
velocities shows a complete decoupling in the system at 
infinitely large Vr. The response of the system in the 
bell-shaped region of the phase diagram of Fig. [2ji leads 
to similar results and can be obtained numerically. Note 
that the Melnikov approximation breaks down for large 
rack velocities and we need a complementary approach 
to examine the behavior of the system in that limit (see 
below) . 

The force-velocity response of the system in the up- 
independent region can also be extracted from this re- 
sult. Figure shows the pinion velocity as a function 
the external load for two values of the rack velocity. 
For Vr < —Vs, the pinion velocity is independent of 
the load forV < |\/P7d - (Vr/R 2 , and drops dras- 
tically after the onset of skipping, whereas the decrease 
upon introduction of load starts from the beginning when 
Vr > ~Vg. One can identify the load at which the pinion 
velocity vanishes as the stall force. 

Strong dissipation. The Melnikov method fails 



if (Vr/R 2 + W 
can proceed by 
v = u in Eq. [TJ 
equation in closed 
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in this limit we 
acceleration term 
allows us to solve the 
which yields u(t) = 
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where r 
The explicit 
to calculate 
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which 



produces very similar curves to those plotted in Fig. 
[5^ above the skipping transition (using the Melnikov 
method). The stall force (load) in this limit can be 



found as W K — F 



1 + (CVr/FR 2 ) 2 - ((Vr/FR 2 ) 



To make a stronger link to the experiments, we need 
to quantify the amplitude of the lateral Casimir force F. 
For a cylinder of radius R located at a (nearest) distance 
H from a plate (see Fig. [TJs), the Casimir interaction 
can be calculated from the corresponding interaction be- 
tween two parallel plates, using the Proximity Force Ap- 
proximation (PEA) 19]. For the normal Casimir force 
between perfect metals, it has been recently shown that 
this approximation works surprisingly well for H < R 
[20| . Using PEA for the lateral Casimir force between a 
pinion of length L and corrugation amplitude a\ and a 
rack of corrugation amplitude with a\,a2 <C H, we 



4 



V s (/Jm/s 




H/X 



R = 200nm, A = lOOnm 
R = lOOOnm, A = lOOOnm 
R = lOOOnm, A = lOOnm 



60 80 100 120 140 160 180 200 



H (nm) 



FIG. 3: Skipping velocity as a function of the gap size for per- 
fect metallic boundaries, corresponding to p = 19.3 gr/cm 3 
(gold) and a\ = a,2 = 10 nm, for different values of radius and 
corrugation wavelength. Inset: the parameter a introduced 
in the text as a function of H/X for perfect metals (red dots) 
compared to the pairwise summation approximation (green 
solid line). The dashed-line shows the empirical formula pre- 
sented in the text. 



the assumption of perfect metallic boundaries, and one 
expects corrections for gap sizes smaller than the plasma 
wavelength of the metals [3] • 

The value of the friction coefficient £ is also instru- 
mental in determining the behavior of the system. While 
this quantity is highly system-dependent in general, it 
is interesting to note that a contribution to dissipation 
also comes from the interplay between the electromag- 
netic fluctuations and the dielectric loss properties of the 
two objects [2lj |. 

In conclusion, we have proposed a design for a nano- 
scale rack and pinion without contact by employing 
the quantum fluctuations of the electromagnetic field. 
This design, and the corresponding variants that could 
be readily conceived, might help towards making more 
durable machine parts for small mechanical system. 

This work was supported by EPSRC under Grant 
EP/E024076/1 (R.G.). 
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to the leading order, where a = t5 ^*_ 1 J (x*) 

with J being the ( "Josephson-" ) coupling function pre- 
sented in Ref. [12|. The inset of Fig. [3] shows the de- 
pendence of a on H/ A for perfect metals and a com- 
parison to pairwise summation approximation, with the 
dashed-line showing an empirical approximate formula 
a e = 1/ cosh 4 / 9 which could be useful for prac- 

tical purposes. Using this result we find 

VS= [-Son ) \W*) \W> [r) 

(5) 

where p is the mass density of the pinion. This result 
shows a strong power law behavior at small values of 
H followed by an exponential decay at large H whose 
length scale is set by A. Figure [3] shows the skipping 
velocity as a function of the gap size, for different val- 
ues of the radius and wavelength. The typical values for 
Vs , which correspond to velocities that the system could 
transfer robustly, are remarkably high: translated into 
angular velocity they are in the kHz region. The strong 
dependence of the skipping velocity on the gap size can 
be used to explore the parameter space and change the 
behavior of the system. For example, one can reduce 
Vs — move in the vertical direction in the phase diagram 
of Fig. [2}j — by increasing H, and thus switch the system 
from the locked-in phase to a reverse gear, at constant 
rack velocity. In other words, changing the separation 
provides a continuous analogue of the clutch gear sys- 
tem. The calculations presented here have been based on 
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